Tunneling phase gate for neutral atoms in a double-well lattice 
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We propose a new two-qubit phase gate for ultra-cold atoms confined in an experimentally realized 
tilted double-well optical lattice [Sebby-Strabley, et ai, Phys. Rev. A 73 033605 (2006)]. Such a 
lattice is capable of confining pairs of atoms in a two-dimensional array of double-well potentials 
where control can be exercised over the barrier height and the energy difference of the minima of 
the two wells (known as the "tilt"). The four lowest single-particle motional states consist of two 
pairs of motional states in which each pair is localized on one side of the central barrier, allowing for 
two atoms confined in such a lattice to be spatially separated qubits. We present a time-dependent 
scheme to manipulate the tilt to induce tunneling oscillations which produce a collisional phase gate. 
Numerical simulations demonstrate that this gate can be performed with high fidelity. 
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Quantum information processing with neutral atoms 
in optical lattices holds great promise, but achieving 
the necessary single-qubit and two-qubit control remains 
an elusive experimental challenge. Recent experimen- 
tal work has made great progress on this front by devel- 
oping techniques to isolate and manipulate an ensemble 
of pairs of atoms in a double- well lattice [j], 0. Previ- 
ous theoretical proposals for quantum gates include di- 
rect collisional interactions in state-dependent potentials 
[H, 0] ! long-range dipole-dipole interactions [H , and Ry- 
dberg states for maximal speed [fjj]. Each involves using 
the internal states of the atom as qubits. 

An alternative is to use the center-of-mass degree of 
freedom, which can be manipulated by either using vi- 
brationally excited states 0, E| , or lattices with less-than- 
unit filling [i| . These proposals require careful control of 
the qubit states and operate in the limit of weak tunnel- 
ing. However, in a recent elegant experiment [lj|, both 
the spatial and internal degrees of freedom of two atoms 
were used to demonstrate the elements of a two-qubit ex- 
change gate. This experiment was in the regime of strong 
tunneling, bringing both atoms into the same well. Tun- 
neling of single and pairs of atoms has been explored in 
an independent double- well experiment (llT | . 

In this paper, we show how a very simple quantum 
logic gate can be implemented using only the vibrational 
state s I12L 1 1 3 1 of the recently realized double- well lattice 
[H, 0, lldj . and using simple control protocols well suited 
to experimental realization. Central to this scheme is 
controlling the double wells to switch on the tunneling 
of an atom from one well to the other. This tunneling, 
in turn, introduces collisions leading to an overall quan- 
tum phase. Furthermore, by using the vibrational states 
this gate does not depend on the internal state, and thus 
requires no active stabilization of the magnetic field. 

A two-dimensional array of double wells is created 
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FIG. I: (a) Two-dimensional optical lattice setup. An ad- 
ditional optical lattice in the z-direction (out of the plane) 
provides three-dimensional confinement, (b) The effective 
potential U(x) (see text) for the double- well lattice with 
e z = tt/2 (black) and Z = tt/2 + 0.5 (gray) for Z f = 0.1 1 
and Vb = 40Er. The potential is plotted in units of the re- 
coil energy Er and as a function of the dimensionless position 
kx. (c) The effective potential U (x) for the double- well lattice 
with Zf = O.f f (black) and Zf = 0.2 (gray) for 6 Z = tt/2 and 
V = 40E R . 



[iL I3| by intersecting at right angles two pairs of coun- 
terpropagating laser beams, with an independent optical 
lattice along the orthogonal direction. This is used to 
trap and manipulate cold 87 Rb atoms. The four beams 
are obtained by folding and retro-reflecting a single laser 
beam as shown in Figflja). All four beams lie in a hori- 
zontal plane. The electro-optic modulator labeled EOM 
/? is used to rotate the polarization of the incoming beam 
[solid line in Fig.QJa)]. Consequently, the sum of the elec- 
tric fields of the four beams in the lattice region has com- 
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ponents in both the horizontal (xy, in-plane) and vertical 
(z,out-of-plane) directions. The EOM devices labeled 8 
and 4> serve to introduce phase shifts (9 Z and 4> z ) between 
these two electric-field components by altering the opti- 
cal path lengths. 

The potential experienced by the atom is proportional 
to the total intensity of the light. In this case the full 
optical potential is 

V(x,y,z) = --Vi[4 + 2cos(2ky + 2<l> xy ) 
+2 cos (2kx — 29 xy — 2(f> xy ) ] 
-V 2 [cos(kx -9 Z - 4> z ) + cos(fcy + <f> z )] 
+V 3 sin 2 (kz), (1) 

where V± and V 2 are proportional to the fractions of 
the total light intensity contributed by the in-plane and 
out-of-plane polarized light, respectively, V3 to the light 
intensity for the optical lattice along the z-axis, and 
k = 2tt/X, where A is the wavelength of the incident 
lasers. It will be convenient in what follows to define 
V and Z f such that Vi =V (1- Z f ) and V 2 = V Z f . 
The electro-optic elements enable control of the relative 
phases 6 Z — 9 xy and <j) z — 4> xy . In this paper we shall only 
consider the case where 9 xy = 4> xy = <j) z = 0. 

When the depth Vq of the lattice is sufficient, the sys- 
tem can be arranged so that each double-well lattice 
site contains exactly two atoms. In describing a single 
two-atom pair, we may consider a one-dimensional "cut" 
through the potential at z — y = and along the x axis, 
U (x) = V{x, 0, 0). Using the three parameters appearing 
in this potential, 9 Z , Zf, and Vq, control can be exercised 
over the "tilt", the barrier height between the two wells, 
and the overall depth of the lattice. The tilt, defined 
as the energy difference between the left and right well 
minima, is controlled mainly by changing 9 Z as can be 
seen in Fig. |Tfb). Tuning the value of Zf changes the 
height of the central barrier, as seen in Fig. QJc). The 
phase gate operation proposed in this paper is performed 
by changing the value of 9 Z only. 

For a single atom in one site of the optical poten- 
tial described by Eq. ([T]), we encounter a number of 
energy scales. By expanding the potential V(x, y, z) in 
y and z we find V(x,y, z) s» U(x) + (\/2)muj 2 y 2 [l + 
Zf cos(kx — 9 Z )] + (l/2)mu^ 2 , with hw y — ^/AErVq and 
fi^z = v^ErVz, where we have introduced the recoil en- 
ergy Er = (hk) 2 /(2m) with m the atomic mass. For 
typical lattice parameters Zf = 0.11, V3 = Vo = 4QEr, 
we find that ftw y — fku z — 12.65-Er. These are the ap- 
proximate energies for excitation along the y or z direc- 
tions, and are larger than an excitation in the double- 
well potential U(x), which is approximately IOEr. In- 
terestingly, the residual coupling between x and y does 
not lead to any first-order mixing of the states, but does 
shift the central barrier by Zf^JVo/En/i w 0.2Er, which 
is small compared to the typical barrier height of 20Er 



(see Fig. [T]) , and thus can be neglected. In the double- 
well potential, pairs of nearly degenerate levels (right and 
left) are actually split. For the ground doublet this split- 
ting is very small, but for the upper doublet, this tunnel 
splitting 2J ps Q.3Er (2J/h is the tunneling frequency). 
As our gate involves only manipulations of energy at the 
level of J, we will restrict our attention to one-dimcnsinal 
motion along the double- well axis (x). 

For two atoms, we must also consider the effective one- 
dimensional atom-atom interaction, which arises from 
the three-dimensional interaction (AttH 2 'a s / 'm)5(ri — r^), 
where a s is the s-wave scattering length of the freely scat- 
tering atoms. By integrating the <5-function over y and 
z, the interaction strength gm along x can be calculated 
(more sophisticated methods are described in (HI)- Us- 
ing the ground state harmonic oscillator wavefunctions 
for y and z we find g\o = 2h^juj y uj z a s . Note that this 
value can be conveniently described by the dimension- 
less parameter gm = kgm/En = 8Tr(a s /\)y/Vo/ER, 
which with typical experimental parameters (a s = 5.3 
nm, A = 810 nm, and Vq/Er = 40) leads to gm ~ 1. 

To design the two-qubit gate operated in this optical 
lattice, we now consider the energy-level spectrum of two 
interacting particles in the double-well potential, using 
the following Hamiltonian: 

2 2 

H = + u (a;i) + £ + u ^ + 9106 (xi _ X2) 1 (2) 

where p± and p 2 are the momenta of atoms 1 and 2. 
Each atom is assumed to move in a single period of the 
double-well potential such as that shown in Fig. [TJ 

We propose to use as qubits the single-particle vibra- 
tional states supported by the double-well potential. If 
the double-well potential is sufficiently tilted and has 
a high-enough central barrier, the lowest four single- 
particle motional-state energies will lie below the barrier 
and will divide into pairs of states with their wavefunc- 
tions well-localized on one or the other side of the barrier. 
The key to our proposed two-qubit gate is to initially 
set the potential parameters such that these particular 
single-particle motional states exist and can be individ- 
ually addressed. Gate operation proceeds by varying the 
tilt (via the 9 Z parameter) in such a way as to entangle 
the atoms, finally ending up with the same double- well 
potential. 

With this goal in mind, we now calculate the eigen- 
values of the two-particle Hamiltonian in Figure [5] as a 
function of A9 = 9 Z — ir/2 (proportional to the "tilt") 
with experimentally relevant parameters. The energies 
of the first fourteen states are shown and labeled by 
their dominant composition (for large tilt) in terms of 
the products of the lowest six single-particle states of 
the double-well: \0 R ), \Q L ), \1 R ), \1 L ), |2_), and |2+). 
The states subscripted by L or R are localized primar- 
ily in the left or right wells, while for these parameters 
|2_) and |2 + ) are not. Also shown next to each of Fig. 
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FIG. 2: Energy levels of the double-well potential, for S7 Rb atoms in a double- well lattice as a function of AG = 9 Z — ty/2 with 
Vo = 40Er, Zf = 0.11, and gm = 1. Only those states with vibrational excitation along the direction of the double-well are 
calculated. While each apparent crossing of eigenvalues is avoided, the first-order tunneling splittings have been highlighted by 
the gray line segments in (b) and (c). For each eigenvalue an approximate label ( |0_r0_r) , etc.) indicates the eigenvector, valid 
for large AO. Next to each panel are schematic energy configurations for each eigenvalue that indicate the occupations of the 
single-particle energy levels. The boxed configurations are the two-qubit states of the system. 



[5]Ja)-(c) are schematic configurations showing which of 
the single-particle levels are occupied. The boxed con- 
figurations correspond to |0l0_r), |1_l0,r), \0l^r), and 
lli/lii), the natural two-qubit states. Note that in the 
above, Bose-symmetrization of the two-particle states has 
been suppressed to simplify the labels (e.g., \0l®r) = 
2- 1 / 2 |0 L )® |0 i? .)+2- 1 /2| 0it )®|0 L )). 

There are several things to observe in the eigenvalue 
spectrum. First, for zero tilt (A9 — 0), there is complete 
symmetry between left and right. Furthermore, those 
states with atoms localized in the same well are shifted 
up in energy. For example, in Fig. [^a) states |0z,0l) and 
|0_r0_r) have greater energy than \0l0r) at A9 = 0. This 
difference in energy is the collisional interaction energy 
Uqo and is approximately one recoil energy (Er) for these 
parameters. Second, for sufficiently large tilt (AO as 0.1), 
the qubit states are both spatially separated and sepa- 
rated in energy, so that single-qubit manipulation can be 
performed on each of the two atoms. 

Between these two extremes, at intermediate tilt, we 
observe several apparent intersections of the eigenvalue 
curves in Fig. [2j Each of these is actually an avoided 
crossing, most due to second-order tunneling_processes 
(J 2 /U) as studied in a recent experiment Uj. These 
splittings are very small and will not be used in the fol- 
lowing. Broader first-order tunnel splittings arc found 
between states |1z,0,r) and |0.r1.r) in Fig. [2jb) and be- 
tween states |1_l1_r) and \1r1r) in Fig. Etc). The corre- 
sponding energy levels have been highlighted in gray for 
emphasis. Note that these two avoided crossings occur at 
slightly different tilts: the first occurs at A8 as 0.05, the 
second at A9 as 0.04. They also have slightly different 
minimum splittings, which we now consider. 

In a simple two-state model of each avoided cross- 
ing, the two splittings are approximately given by 



y/(U 01 - AE) 2 + 4J 2 and y/(Un — AE) 2 + 8J 2 , where 
2 J is the single-particle tunnel splitting (of the upper 
doublet; the second splitting is enhanced by a factor of 
y/2 due to Bose symmetry), AE is the energy tilt, and 
Uqi and U\\ are the interaction energies for the two states 
|0_r1_r) and \1r1r), respectively. Note that, to achieve 
maximal tunneling of the atoms, one must tilt the dou- 
ble wells to AE ~ U, using the tilt to compensate for 
the interaction energy. The two interaction energies are 
different; using a harmonic oscillator approximation, we 
find Uoi = Uoo, while Un — 3E/oo/4. However, by choos- 
ing the tilt such that AE as 7£/"oo/8 — 8 J 2 /Uoo, one can 
induce tunneling (from left to right) with the same oscil- 
lation frequency starting from either of the qubit states 
or|lilfl). Both interactions and strong tunneling 
are required to achieve this type of sychronized oscilla- 
tion. 

We now turn to the gate sequence, illustrated in FigJ3] 
Starting from large tilt (time 1), in which the two-qubit 
states can be individually addressed, we shift Ad to an 
intermediate tilt optimized such the tunneling frequen- 
cies are equal. Holding the lattice at this tilt for some 
amount of time (including time 2), collisions will occur 
when both atoms occupy the right well. After some time 
an appreciable phase will develop, after which a tilt back 
will recover the original two-qubit states (time 3), with 
an overall controlled phase. 

The synchronized tunneling oscillations are demon- 
strated in FigJ3Jb), by solving the time-dependent 
Schrodingcr equation with the Hamiltonian ([2]) using a 
split-operator spectral method and system parameters 
from Fig. [2] Here the atom-atom interaction has been 
approximated by a superposition of two Gaussians, opti- 
mized to match the energy levels shown in Fig. [5J Only 
the tilt A8 — 8 Z — tt/2 is manipulated, using a ramp from 
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A8i = 100 mrad to A8h — 33 mrad in time t r = 0.12 ms, 
and a hold time of th — 1-46 ms, for an overall gate 
time of t g = th + 2t r = 1.7 ms. The two-particle prob- 
ability density for various steps during this sequence is 
shown in FigJ3{c). To interpret these panels, observe 
that the initial state has little amplitude near the diag- 
onal (x\ — X2), indicating that the atoms are well sepa- 
rated. During the hold, when one atom tunnels from left 
to right, there is a large probability along the diagonal, 
indicating that both atoms are in the right well. Finally, 
at the end, the wavefunction has returned to (approxi- 
mately) the initial condition. 

After two oscillations of the probability density, the 
overall controlled phase that is accumulated (not shown) 
is (ft = <f)n + O o — 0oi — 0io ~ 0.9-7T (this phase can be 
adjusted by changing the ramp parameters; see below). 
To quantify the accuracy of this quantum logic operation, 
we construct the two-qubit gate's matrix representation 
V by evolving each of the initial states corresponding to 
\0l0r), |1l0^), |0l1_r), and |lz,l fl ), and after removing 
single- qubit phases we calculate the average gate fidelity 
T til: 

^^EK^^HM 2 '^) (3) 

j,k=0 

where W = diag(l, 1, 1, e l< ^) and cn are the Pauli ma- 
trices for the first qubit. This is equivalent to an average 
of the gate fidelity over all initial two-qubit states. 

This fidelity is calculated for various values of the ramp 
and hold times t r and th in Fig. H](a), and for various val- 
ues of the final tilt A9h and hold time in Fig. [U(b). Even 
for this simplest of control sequences, fidelities greater 
than 0.99 are possible. Furthermore, by operating in the 
limit of strong tunneling, controlled by the tilt, we have 
achieved a gate operation that is nearly ten times faster 
than previous proposals [3, HI]. While promising, this gate 
is not quite optimal: a gate equivalent to a controlled- 
NOT requires <p = n. Improvement of both the controlled 
phase and overall fidelity should be possible by waiting 
longer, using tighter transverse confinement, a deeper lat- 
tice, or through optimal control techniques [T3|. Other 
issues that may limit the fidelity are coupling to vibra- 
tional excitations in the transverse directions, and will be 
studied elsewhere. Nevertheless, Fig. 2] shows that small 
imperfections in the control parameters (due to timing 
errors, or inhomogeneities across the lattice) still allow 
for fidelities greater than 0.9, very promising for initial 
experimental demonstration. 

In summary, we have analyzed the interaction of two 
neutral atoms in one cell of the double-well optical lat- 
tice. By manipulating just one property of the potential, 
the "tilt", we have shown how to achieve a high-fidelity 
controlled phase gate. By operating in the regime of 
strong tunneling, this scheme leads to a fast gate oper- 
ation and is experimentally accessible without the need 
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FIG. 3: (a) Schematic implementation of the tunneling phase 
gate consists of a three step sequence, from large (1) to small 
tilt (2) and back (3). During the hold period (including time 
2), an atom in the first excited state of the left well will tunnel 
to and from the right well, (b) Time oscillations of the qubit 
populations given the initial condition j*(0)) = \1l0r) (solid) 
and |\P(0)) = 1 1 i 1«) (dashed). Also illustrated (gray solid) 
is the optical lattice parameter A9 = 9 Z — n/2 (right axis), 
(c) Two-atom probability densities |\P(a;i, 352) | 2 , for various 
stages of the tunneling phase gate sequence, as indicated in 
(b), and the initial condition |\P(0)) = |1l1.r}. 



to control the internal states of the atom. Consequently, 
this protocol can be performed without stabilizing the 
magnetic field environment, avoiding a significant source 
of decoherence. 
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The fidelity maximum is 0.994. (b) Average gate fidelity T 
as a function of the hold time t h and the final tilt A9f, with 
t r = 0.12 ms, A9i = 100 mrad, V = 40E R , and Z f = 0.11. 
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